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Abstract 

We analyze a perturbation of the boundary Sine-Gordon model where two boundary 
terms of different periodicities and scaling dimensions are coupled to a Kondo-like 
spin degree of freedom. We show that, by pertinently engineering the coupling with 
the spin degree of freedom, a competition between the two boundary interactions 
may be induced, and that this gives rise to nonpertubative phenomena, such as the 
emergence of novel quantum phases: indeed, we demonstrate that the strongly cou- 
pled fixed point may become unstable as a result of the "deconfinement" of a new 
set of phase-slip operators -the short instantons- associated with the less relevant 
boundary operator. We point out that a Josephson junction network with a perti- 
nent impurity located at its center provides a physical realization of this boundary 
double Sine-Gordon model. For this Josephson junction network, we prove that 
the competition between the two boundary interactions stabilizes a robust finite 
coupling fixed point and, at a pertinent scale, allows for the onset of 4e supercon- 
ductivity. 
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1 Introduction 



There is a large number of pliysical systems tliat can be mapped onto quantum 
impurity models in one dimension [1]. Embedding a quantum impurity in a 
condensed matter system may alter its responses to external perturbations [2] , 
and/or induce the emergence of non Fermi liquid, strongly correlated phases 
In quantum devices with tunable parameters impurities may be realized 
by means of point contacts, of constrictions, or by the crossing of quantum 
wires or Josephson junction chains [HI5|5|7] . While a standard perturbative 
approach works fine when impurities are weakly coupled to the other modes 
of the system (the "environment"), there are situations in which the impuri- 
ties are strongly coupled to the environment, affecting its behavior through 
a change of boundary conditions: when this happens, it is impossible to dis- 
entangle the impurity from the rest of the system, the perturbative approach 
breaks down, and, consequently, one has to resort to nonperturbative meth- 
ods, to study the system and the impurity as a whole. Such nonperturbative 
tools are naturally provided by boundary field theories (BFT) [T||8] : BFTs al- 
low for deriving exact, nonperturbative informations from simple, prototypical 
models which, in many instances, provide an accurate description of experi- 
ments on realistic low dimensional systems [9]. In particular, BFTs have been 
successfully used to describe Josephson current patterns in Josephson devices, 
such as chains with a weak link [TOpTj . SQUIDs [T2IT3] and Y junctions [7]. 

Motivated by the Kondo effect [H], impurity models have been largely stud- 
ied to describe some magnetic chains [15] , and static impurities in Tomonaga- 
Luttinger liquids (TLL)s [16]. A renormalization group approach to those sys- 
tems leads, after bosonization [T7], to the investigation of the phases accessible 
to pertinent boundary sine-Gordon models [16j. Scattering from an impurity 
often leads the boundary coupling strength to scale to the strongly coupled 
fixed point (SFP), which is rather simple since it describes a fully screened 
spin in the Kondo system or a severed chain in the Kane-Fisher model [18]. A 
remarkable exception is provided by the fixed point attained in overscreened 
Kondo problems, where an attractive finite coupling fixed point (FFP) emerges 
in the phase diagram [13] ; this FFP is usually characterized by novel nontriv- 
ial universal indices and by specific symmetries. In the analysis of the Kondo 
effect, an SU{2) invariant coupling of a local spin degree of freedom with the 
spin density of conduction electrons, allows for engineering a marginally rel- 
evant interaction, which would otherwise be irrelevant. Similar behaviors are 
realized with crossed TLLs where, as a result of the crossing, some operators 
turn from irrelevant to marginal, leading to correlation functions exhibiting 
power-law decays with nonuniversal exponents [T91I6] . 




Superconducting Josephson devices allow to engineer remarkable realizations 
of the above situations, [TTP^] . For superconducting Josephson chains with 
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an impurity in the middle p^p!T] or for SQUID devices [T2|T3] the phase di- 
agram admits only two fixed points: an unstable weakly coupled fixed point 
(WFP), and a stable one at strong coupling, while, for pertinent values of the 
fabrication and control parameters, a FFP emerges in Y-shaped Josephson 
junction networks (JJN)s |7]. The boundary field theory approach developed 
in Ref. |ll|IT3] not only allows for an accurate determination of the phases 
accessible to a superconducting device, but also for a field-theoretical treat- 
ment of the phase slips (instantons), describing quantum tunneling between 
degenerate ground-states; furthermore, it helps to evidence remarkable analo- 
gies with models of quantum Brownian motion on frustrated planar lattices 

EQEU. 



Here we study the effect of adding a less relevant scaling operator to a bound- 
ary Sine-Gordon model. Most analytical computations hold only when the 
second less relevant operator has been scaled away [22]: conventional wisdom 
suggests indeed that one should be able to neglect all less relevant operators, 
when computing properties close to the infrared fixed point. However, this 
expectation is based only on weak coupling expansion, which can be quite 
misleading [23|2l] . In this paper, we shall exhibit an explicit example of a 
boundary field theory model where the added perturbation may become rele- 
vant at strong coupling and we shall provide a superconducting device where 
the onset of new nonperturbative phenomena may be observed. Adding to a 
boundary Sine-Gordon model a perturbation with a different scaling dimension 
and periodicity allows, in a superconducting device, to change the tunneling 
charge and, thus, to affect the transport across the device. For quantum Hall 
fluids [25], superconductor- normal metal contacts [25] and Kondo quantum 
dots [2Z], adding a perturbation modifles the charge of the excitations , as 
evidenced in dc shot noise measurements [281I29II30I. 



We shall consider a boundary fleld theory with two boundary terms, of differ- 
ent periodicities and scaling dimensions, coupled to a Kondo-like spin degree of 
freedom. The resulting model is described by a boundary double Sine-Gordon 
(BDSG) Hamiltonian, given by H^bsg = -f^LL + Hb, where i^LL is a spinless 
one-dimensional Tomonaga Luttinger Hamiltonian [21] -deflned on a support 
of length L, with velocity u and Luttinger parameter g- given by 



L 





1 ^ 
u [~dt 




and 

Hb = -9i cos[$(0)] - g2 cos[2$(0)] - 5|,S" - 5^S^' , (2) 
describing the interaction between the Luttinger fleld $ and a spin- 1/2 de- 
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gree of freedom, localized at a; = 0. In this paper, we shall show that one can 
engineer the coupling with the spin degree of freedom, so as to induce a com- 
petition between the two periodicities in H^, leading, in some instances, to the 
emergence of new quantum phases. We shall show indeed that, for 1 < g < 4 
and for B\\ = B± = 0, the less relevant interaction —(72 cos [2$ (0)] destabilizes 
the strongly coupled fixed point, as a result of the "deconfinement" of new 
phase-slip operators (instantons), characteristic of the double Sine-Gordon in- 
teraction [32j. To fix the ideas, we analyze in detail the Josephson junction 
network depicted in FiglH since it provides a remarkable physical realization 
of the BDSG model described by H^; in this JJN, we show that the competi- 
tion between the two periodicities in stabilizes a robust [33fM] FFP, and 
-at a pertinent scale- allows for the emergence of 4e superconductivity [35] . 

The paper is organized as follows: 

In section [21 we show that the JJN in Fig{T]is indeed described by H^, i.e., by 
a double boundary Sine-Gordon Hamiltonian coupled to a pertinent spin-1/2 
local spin degree of freedom; 

In section [3l we determine the phase diagram of the DBSG model, using the 
renormalization group (RG) approach and show that it admits a WFP, a 
strongly coupled fixed point (SFP), and, for 1 < < 4 and for = B± = , 
a FFP. Furthermore, we show that, near by the FFP, the emerging local spin 
degree of freedom is robust against decoherence; 

Section mis devoted to the analysis of Josephson current patterns exhibited by 
the JJN. There we show that 4e superconducting correlations may be probed 
in a Josephson current measurement, in all the phases accessible to the JJN; 

In section [5l we evidence that a shot noise measurement can account for the 
emergence of 4e tunneling charges in the JJN, near by the WFP. Furthermore, 
to show that 4e superconductivity is a feature of the JJN also far from the 
WFP, we derive an exact formula for the dc current, as well as for the shot 
noise, at the "magic point" g = 2 [36j, where the WFP is not IR stable; 

Section [H] is devoted to our concluding remarks, while the appendices provide 
the necessary mathematical background for the analysis carried in the paper. 



2 The boundary double Sine-Gordon Hamiltonian 

In this section, we show that the JJN depicted in Fig{Tl may be effectively 
described by Hbdsg, defined in Eqs. (fTp|) . In Eq.([2]), gi and g2 are real param- 
eters, with g2 > 0, while and are, respectively, the z-component and the 
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Fig. 1. The Josephson device: the central region C with two leads connected to two 
external contacts. 

x-component of a spin- 1/2 operator. and B± may be regarded as the two 
components -along z and x, respectively- of an external magnetic field acting 
on S; as such they may be regarded as control parameters to tune the onset 
of different regimes. 

The spin- 1/2 degree of freedom allows for i^BDSG to be invariant under 



n : 



(3) 



which realizes the usual "Sine-Gordon symmetry" with period vr (instead of 
2tt); for k odd, ri involves also the sign inversion of 0. As we shall see, the 
emergence of this symmetry is crucial to account for the novel behaviors in 
the JJN depicted in FigHJ 

The JJN consists of a central rhombus C, made with four Josephson junctions 
of nominal strength J, pierced by a dimensionless flux (p (i.e., (p = $/$o) 
connected to two chains (leads) of Josephson junctions, of nominal strength 
Ej, with charging energy Ec ^ Ej, and charge repulsion strength between 
nearest-neighboring junctions given by E^. The gate voltage applied to each 
junction is tuned at the degeneracy between charge eigenstates with N and 
N + 1 Cooper pairs, so that each junction may be regarded as an effective 
spin-1/2 variable. In this regime, the low-energy, long wavelength dynamics 
of the two leads is well described in terms of two LL Hamiltonians for the 
plasmon fields of the chain on the left- and the right-hand side respectively, 
$<,$>: the Luttinger parameters q and u are given by gf = — — -7^, 

ox- ^ o J ^ 2(7r-arccos(-Y)) ' 



^ Notice that this is consistent with keeping unchanged, as one may change sign 
to two components of S, say S^, (not appearing in H^), without altering the 
canonical commutation relations 
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U = Vf 



2 arccos(^) 



(A = {E' - 3Ej/16E^)/Ej) [TOlITT]. The central region C 

is described by He = - J E?=o{e'^5+57+i + h.c.} + J" E?=o with 5,- 

spin-1/2 variables defined at site j . To trade He for an effective boundary 
interaction, one performs a systematic Schrieffer- Wolff (SW) sum over the 
high-energy eigenstates of He- This is carried out in appendix lAl where it is 
shown that, for ip = n, the ground state of He is twofold degenerate, with 
the two degenerate states given by | ft) = i^llv^ll titi) + I Itlt)] + I iitt 

) + 1 tm) + 1 mt) + 1 ;m)}, and by u) = ^{v^[i nu) - 1 mt 
)] - ^1 un) - ^1 tm) + ^1 mt) + tttt)}- 

The degeneracy between | ft) and | Jj.) may be removed by slightly detuning 
(p to (p = IT + S, with |5/7r| ^ 1. From Eq. flA.20p . one sees that removing the 



degeneracy induces the term — i?|jS^ in H^; at variance, detuning the gate 
voltage applied to the junctions in C yields the term —B±S^ (see Eg. flA.2ip ). 

Connecting the leads to C with two Josephson junctions of nominal strength 
A Ej, J), allows -via the SW procedure described in appendix |A]- to de- 
termine the effective boundary Hamiltonian, which, to the fourth order in 
A, coincides with Eq.®, with = -4Jsin(^), B± = ^^^z^^^, 
$ = $>)/V2, gi = § (ir^)' 92 = 2^, C being a numerical co- 



efficient ~ 10^ . As evindenced in Ref . |35] . the first term in H^ describes 
tunneling of Cooper pairs between the two leads of the device, while the sec- 
ond term is responsible for the coherent tunneling of pairs of Cooper pairs 
across C. 

While -f^B provides the dynamical boundary conditions (BC)s at the inner 
boundary, the BCs at the outer boundary (x = L) depend on the type of ex- 
ternal contacts one attaches to the JJN to induce a current across the leads. In 
particular, when the JJN is connected to two metallic leads at a finite voltage 
bias V (or not contacted), one may safely assume Neumann BCs (^|^ = 0) 
at the outer boundary while, when the device is connected to two bulk su- 
perconductors at fixed phase difference a (as it happens when a dc Josephson 
current is induced across C), one may assume Dirichlet-like BCs {^{L, t) = a) 
at the outer boundary. 



3 Perturbative renormalization group analysis 



In this section, we use the RG approach to investigate the phase diagram 
accessible to a system described by i^DBSG- A perturbative analysis of the 
boundary interaction shows that there is a range of values of g for which H^ 
becomes a relevant operator and, furthermore, evidences the effects of its two 
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competing harmonics. We shall show that the phase diagram admits a WFP, 
an SFP, and, for 1 < (7 < 4 and for = = 0, it allows for the emergence 
of a FFP, which is responsible for some of the remarkable novel behaviors 
exhibited by the JJN in FigJTJ 

In subsection 13.11 we derive the perturbative RG equations at weak coupling 
and use them to investigate the stability of the WFP; in subsection 13. 2[ we 
repeat the same analysis near by the SFP and, in subsection 13. 3[ we show 
that a stable finite coupling fixed point emerges within a pertinent window of 
values of g when B\\ = Bj_ = 0; finally, we investigate how decoherence may 
be frustrated [53] when the JJN is operated near by the FFP. 



3.1 Perturbative renormalization group analysis near by the WFP 



To derive the perturbative RG equations near by the WFP requires computing 
the partition function of the system by integrating over the field $, as well as 
over the spin variable S. To perform the latter integration, we resort to the 
imaginary time formalism and introduce two local complex fermion variables, 
{a, b}, to describe the spin-1/2 operator, which is then given by — > a^a — 
b^b, and — > a^b + b^a. As a result, the Euclidean action for the boundary 
degrees of freedom of the system is given by 

S^ = SS^ + S^^ , (4) 



where 



0(0) 



dr < 



— - za;o - 5,1 



+ b^ 



9 ■ 

— -iuj^ + B\\ 



b\-B±l dr {a^b + b^a} , (5) 



and 



-gi j (ir {a'l'a - cos[$(r)] - 5(2 j dr cos[2$(r)] 



(6) 



with Wo = vr//3, /3 = {kBT)-\ and $(r) 
function is then given by 



Z = Zo (T,e-^B 



$(0,t = ir) The partition 



(7) 



with Zq = Tr exp[— /3(ifLL ~ ^W^^ ~ B±S^)], being the imaginary time- 
ordering product operator, and (...) denotes thermal averaging with weight 
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function exp[-/3(ffLL - SyS^ - BxS^)]/Zo. 

To integrate over the local fermion operators, one needs to determine the 
relevant imaginary time correlation functions of S^, S^; these are given by 

(S^)(o) = cos(^) , (S^)(o)=sin(^) , (8) 



with cos(^) = B./JBl + Bl, sm{9) = B^/ JBf^ + Bl, and 



(T. [S^ (r) (r' )] ) (0) = cos^ (6) + sin' (6)6-''^^-^' I 
(T, [S" (r) (r )] ) (o) = sin^ (6) + cos^ {e)e-'^\^-^' I 



(9) 



with A = + Bl. From Egs.®, one sees that 5j| = (5^ = 0) ^ cos(0) = 
0(sin(^^)=0)^(S^)(o) = 0((S-)(o) = 0). 

To integrate over $, one has to specify its BCs at both boundaries. At x = 0, 
the pertinent BCs are set by energy conservation, which amounts to require 



27r dx 



(7iS^sin[$(r)]+2(?2sin[2$(r)] . 



(10) 



Within a perturbative approach in H^, one should then assume Neumann BCs 
at X = (i.e. 



' dx 



0), and require free BCs at x = L (i.e., ^^1^'""-* = 0). 



Neumann BCs at both boundaries yield the following mode expansion for 
$(x, r) 



$(a;,r) = ^ 



2muT ~ 
(l)o + —^P 



+ ^ E cos 
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n 



2Tmx 



e L 



with [00, -P] = i, and = n6n+m,o- By substituting Eq.l ITT]) into Eq.(j6]), 

and normal-ordering the vertex operators with respect to the ground state of 
i^LL, one gets 



J dr {-g^[a\r)a{r) - h\r)h{r)] : cos[$(r)] : +^2 : cos[2$(r)] :} , (12) 



with gi = (^j^^ " gi, and g2 = ) ® g2, cl/u being a pertinent short (imaginary 
time) distance cutoff. 
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From Eq. (fT2|) . using the standard factorization formula of vertex operators 
[37] . one gets that 



(T,[: e^"i*("^) : . . . : e*""*^"'^) :])o = exp 



E 2^7.(r.,r,) 
i<j=i 9 



with 



7t 



In 



2 sinh -p-(T — r 



(14) 



Eq.(|2]), imphes that, in computing Z as a power series in , vertex operators 
: e^**^"^^ : should be always accompained by an operator S^(r). As a result, as 
long as cos(6') ^ and A 7^ 0, and for length scales L > Lx one finds 

that 

13 

Z = Zo E A / ^^.] E n ( ^) (T.[: e-*(-) : . . . : e-"*(-) :])(o) 

n=0 j=l ^ a,={±l,±2}^=l ^ ^ ^ 



with 



1/1 = cos(6')5(i , y2 = g2 + sin (6*) 



2A 



nu 
4^ 



a \ 9 9 



(16) 



Eq. f|T6l) allows to infer the RG flow of the boundary interaction, when (cos(^), sin(^^)) 7^ 
(0, 0). It should be noticed that Z in Eq. flT^ may be regarded as the partition 
function of a one-dimensional Coulonb and that the fugacities yi, 2/2 associated 

_i 

to its "charges" scale as L 9 and as L 9 , respectively. 

To study the behavior of the boundary interactions along the RG trajectories, 
one needs to write down the RG equations for the running coupling strengths 
Gi{L) = gi (^j^^, G2{L) = 92 where Lq is a reference length; these 

equations may be derived from the short (imaginary time) distance operator 
product expansions (O.P.E.)s of the vertex operators entering Hb [H], which 
are given by 



{: e'^**^'^^ : : e"^**^^ ^ :} ^ _ 



{:e 



7TU[T — T ) 

L 

tcu{t — t') 
L 



. g±2i$(r) . _^ _ _ _ 



: e 



=Fj<I>(r) 



:+... , (17) 



where the ellipses denote non diverging contributions. As a result, one gets 
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\ G2 



L/L 







L/L 







Fig. 2. Renormalization group trajectories of Gi{L) and G2{L), for 1 < (7 < 4 
and for = i?_L = 0, for different values of Gi(Lo) and G2{Lq). Top left panel: 
Gi(Lo) / 0, G2(^o) = 0; Top right panel: Gi(Lo) = G2(^o); Bottom left panel: 
Gi(Lo) < ^2(^0); Bottom right panel: Gi(Lo) = 0,02(^0) / 0. 



dGi 



rfG2 



din 



1 

1 - - 

9, 

4' 

1 - - 
9 



Gi + G1G2 = /?! [Gi, 6*2] 



G2 + - 



1(^1 



2 2 



l^2[Gl, G2] 



From Eq.(lT8l). one sees that, for < 1 both /9i and (32 are negative for small 
values of Gi, 6*2; thus, ioi g < 1, i^B provides an irrelevant boundary pertur- 
bation, and a perturbative analysis in the boundary coupling is expected to 
yield consistent results. On the contrary, for g > 1, Pi starts positive, and, 
despite the fact that (32 starts negative form small values of Gi{Lq),G2{Lq), 



1)P 

'91 



since Gi 



(32 becomes positive at some intermediate scale L12 ^ 
increases as L/Lq increases. Thus, Gi renormalizes G2, so as to make, even 
for g < 4:, cos[2<l>(0)] a relevant perturbation, despite the fact that the linear 
term of the P2 function is negative. 

Upon defining the "healing length" L^, by Gi{L^) ~ 1 IZ1, one finds that, 
for B\\ = Bj_ = 0, and for g > 1 and L > L^,, the system is driven out 
of the perturbative regime, still preserving its fundamental periodicity $ — )■ 
$ + TT (i.e., Ti-symmetry). As we shall see, ri-symmetry allows to classify the 
symmetries of the various infrared (IR) stable fixed points accessible to the 
JJN. 



^ Estimating L^: requires, in general, to resort to a numerical calculation. For a 
simple choice of the bare value of the couplings, such as Gi{Lq) ~ G2{Lq) ~ 2y, 



one obtains 



1 



In 



i-(i/g)+y 



(2-(l/9))!/ 
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3.2 Renormalization group analysis near by the SFP 



In the previous subsection, we showed that, for g > 1, the WFP becomes IR 
unstable and, that, for L > L^,, the system is driven out of the (weakly coupled) 
perturbative regime. To infer the IR properties, one may safely assume that 
the boundary interaction is driven all the way down to the SFP and, then, 
exploit standard boundary conformal field theory (BCFT) techniques [6], to 
build the leading boundary perturbation at the SFP, in order to derive the 
pertinent RG equations. 

When 7^ and 1 < < 4, the running coupling Gi{L) grows with L, while 
G2 is an irrelevant coupling. As a result, for L > L^, the system flows to a SFP 
where either $(0) = 2TTk, or $(0) = 7r(2A; + 1) (A; G Z), according to whether 
cos{6) > 0, or cos(6') < 0. Nothing changes when g > 4, since, even if G2{L) 
becomes a relevant coupling, the ratio Gi{L)/G2{L) diverges, as L — )■ 00. At 
variance, when = 0, the SFP is reached when $(0) = 2Trk, = 1, or 
$(0) = 7r(2fc + 1), = —1 {k integer). In the following, we shall derive 
which is the leading boundary perturbation, at the SFP. 

To construct H^, one starts from the explicit form of ^{x, r) obeying Dirichlet 
BCs both at X = 0, and at x = L, since this BCs allow to determine the zero- 
mode contribution. For this purpose, we need only to set r) = a, since 
$(0) is determined by the boundary interaction. Thus, one gets 



t) = a 

V9 



L — X 
L 



TT 



P 




sm 



TTUX 



(19) 



When 7^ 0, the eigenvalues of the zero-mode operator P are either given 
by y/gi^k + a/vr), or by ^/g{2k + 1 + a/n), according to whether cos(^) > 0, 
or cos(^^) < 0; at variance, when cos(^) = 0, the eigenvalues are given by 
Pk = y/9{k + a/n). Knowledge of the spectrum of P allows for building 
the leading boundary perturbation at the SFP, using the delayed evalua- 
tion of boundary conditions (DEBC) technique, developed in Ref. [B]. Within 
the DEBC approach, a generic boundary perturbation at the SFP may be 
represented as a linear combination of boundary vertex operators Va^b ='■ 
exp[z(a$(0) + 60(0))] :, with Q{x, r) being the dual field of ^{x, r) H], whose 
mode expansion is given by 



Q{x,t) = y/g 



muT 
L 



P 



E- 

^ n 



cos 



nnx 



(20) 



0(x,r) is defined by the cross derivative relations, — q^' ^ 

_ i d1>{x,T) 



de{x,T) 

dx 



i de{x,T) 

u dr ' 



and 
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Due to the Dirichlet BCs at x = 0, $(0) does not enter Va^b, and, from the com- 
mutator [P, : exp[i69(0)] :] = 2b^/g : exp[i60(O)] :, one sees that : exp[ibQ{0)] : 
shifts pk to pk + '^y/gb. Moreover, since pk = y/g{k + a/n), one gets b = m/2. 



m 



6(0) 



m G Z. As a result, in Hb there will be linear combinations of cos 
Since the symmetries of require that a shift in p^ is accompained by a 
change in the sign of S^, odd-m terms must be multiplied by S^. As a result, 
Hb may be written as 



-/iiS^' cos 



9(0) 



/i2 cos[0(O)] -5xS^' -5|jS^ 



(21) 



Remarkably, the boundary interaction Hamiltonian Hb in Eq. ( 12 ip is the same 
as i^B in Eq.([2]), provided that: 



$ o e/2 , ^ S^' , gi,g2 ^ /Ui,/i2 , o 5j 



(22) 



Thus, the RG analysis near by the SFP may carried out by mainly repeating 
what has been already done near by the WFP. 



Normal ordering of Hb is accounted for if one defines /xi = (f ) ^ yUi, = 
(f )^ Furthermore, for B± = 0, repeating the same steps of subsection 13. H 
yields the RG equations for the running coupling strengths Ai(L) = Lfii, and 
A2{L) = Lfl2, which are given by 



dAi 



din 

dA2 



1-^ 

4 



Ai + A1A2 



A2 



(23) 



From Eqs.( |9]l22l) . one sees that, if B\\ ^ 0, the vertex operators : e^2®W : 
are not a relevant perturbation to the SFP. These operators may be regarded 
as the field theory representation of the instanton/antiinstanton trajectories 
between minima of H-qy^sg'-, in comparison with the instanton/antiinstanton 
trajectories associated to the operators : e^*®^*^-* :, we refer to them as "short 
instantons" (SI). Since the correlator (Tt-[: e5®(°) : : e~5®('') is always 

accompained by a correlator of the local spin S^, one sees that, for ^ 
and -B_L = 0, Sis are always confined, and thus associated to a "jump" of 
27r. At variance, for 1 < g < A and B\\ = 0, Sis are deconfined, allowing for 
tunneling between minima separated by vr. Finally, for 1 < g < 4 and for 
B\\ = B± = 0, the system flows towards an attractive, IR stable finite FFP. 
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We plot the phase diagram for B± = B\\ = in Figj3l there we show that the 
boundary interaction flows towards weak coupling for g < 1, towards strong 
coupling for > 4 while, for 1 < < 4, it flows towards a finite coupling 
fixed point. For the sake of completeness, the phase diagrams corresponding 
to different values of B\\ and B± are reported in FigJH From FigJU one sees 
that: for B\\ ^ 0, 5^ 7^ 0, the boundary interaction is relevant or irrelevant 
whether (7 < 1, or (7 > 1 and, since B\\ 7^ always confines the Sis near 
by the SFP, the WFP is IR stable for ^ < 1, while the SFP is IR stable 
only for (7 > 1; at variance, for B^\ = 0, B± 7^ 0, charges ±1 are confined 
over an imaginary time scale ~ A^^ (see Eq.([n])), so that the leading effective 
boundary perturbation at the WFP has scaling dimension 4/g and the Sis are 
deconfined at the SFP. This makes the WFP IR stable for g < 4:, and the SFP 
IR stable for g > 4] for B\\ ^ 0, B^ = 0, one gets the same phase diagram as 
for 5|| ^0,B^y^ 0. 

SI deconfinement is the key mechanism destabilizing the SFP, when the WFP 
is IR repulsive; for this reason, a detailed account of this mechanism is pro- 
vided in appendix [Bl However, a similar - this time perturbative- mechanism 
holds near the WFP since, as long as B\\ ^ 0, charges ±1 are always more 
likely to appear than charges ±2. As a result, while for g < 1 Coulomb gas 
charges are suppressed since the corresponding fugacity is an irrelevant op- 
erator, for 1 < < 4, charges ±1 proliferate, while the fugacity for charges 
±2 still is irrelevant. Eventually, both charges may proliferate for (? > 4, but 
the fugacity for charges ±1 is still more relevant than the fugacity for charges 
±2. Accordingly, for B\\ 7^ and for g > 1, the elementary charged excita- 
tions tunnelling across C carry charge ±1 (in units of the Cooper pair charge 
e*). At variance, for B^\ = 0, the fugacity yi goes to zero and only charges 
±2 go across C, via coherent tunnelling of pairs of Cooper pairs, induced 
by the boundary interaction H^; it should be noticed that, in this range of 
parameters, coherent tunnelling of Cooper pairs near the WFP emerges as a 
perturbative phenomenon. We observe also that, for B\\ = 0, B± 7^ 0, the fu- 
gacity yi goes to zero and only charge ±2 excitations cross C as a result of the 
coherent tunnelling of pairs of Cooper pairs, induced by the boundary inter- 
action i^B- Finally, for B^\ = B± = 0, the fugacity yi is relevant for > 1 and 
the elementary excitation carries charge ±1; as we shall discuss in more detail 
in section m for B\\ = B± = 0, the JJN may be regarded as a " 4e-deconfined 
superconductor" , since it exhibits 2e elementary charged excitations while its 
ground state correlations still support 4e superconductivity. 

In summary, our analysis evidences that a FFP emerges only when the Sis are 
a relevant boundary perturbation to the SFP. As shown in appendix [Bl this 
happens only over length scales L such that L^, < L < L^top ~ Lq i — j . 



For B\\ — )■ 0, Lstop 00, and, thus, the FFP becomes a stable, attractive fixed 
point. 
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Fig. 3. RG flow for = = B± = 0. The boundary interaction flows towards 
weak coupling for g < 1, towards strong coupling for 5 > 4, while, for 1 < g < 4, the 
RG flow points towards a finite coupling fixed point, represented as a gray circle. 
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Fig. 4. RG flow diagram for various values of the Luttinger parameter g and of 
the control parameters B^^,B±: Left panel: RG flow diagram for 7^ 0,B± 7^ 0; 
Middle panel: Flow diagram for = 0, B± 7^ 0; Right panel: Flow diagram for 

Bn / 0,B_L = 0. 



3. 3 Frustration of decoherence near by the finite coupling fixed point 



In this subsection, we focus on the analysis of the interaction between the 
spin-1/2 degree of freedom sitting at x = (S), and the plasmon modes of the 
leads, which may be regarded as a bath of this two-level quantum mechanical 
system. We shall show that, near by the FFP (i.e., for 1 < g < 4 and for 
= 0), the decoherence induced by the plasmon bath is drastically reduced 
[33|M] . To do this, we compute the spectral density of states (SDOS) of the 



local spin-1/2 variable, which is given by 



s{n) 



(24) 
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where is the transverse spin susceptibihty, defined as 

oo 

[X^m = J t/re^^MT.[{cos(^)S^(r)-sin(^)S^(r)}{cos(0)S^(O)-sin(0)S^(O)}]) .(25) 



To carry out our task, it is most convenient to compute for g = 1,4, 

and then to extrapolate the behavior of S{Q) for any 1 < g < A. 

For g = 1, cos[2$(0)] is an irrelevant perturbation and one may safely assume 
that the FFP should lie at a distance 0{gi) from the WFP, so that Neumann 
BCs may be imposed on $ at x = 0. As a result, one gets the simplified 
boundary action given by 

= -^Jdr {a^a - b^b}[: e**(") : + : e'^^^") :] , (26) 





with the operators : 6^**'-°^ : having scaling dimension equal to 1. Similarly, for 
g = A, : e^®*^"-* : has scaling dimension 1, and the Euclidean boundary action 
is given by the dual of Eq. ( I26l) , namely 

= 5^0) -^J dr {a^b + b^a}[: e^®W : + : e'^^^^) :] . (27) 





By means of a Bogoliubov transformation to the normal modes of the bound- 
ary fermions a,/3, one finds 



S'' = -sm{e){a^a-/3^/3} + cos{e){a^P + /3^a} , (28) 

with 

a = cos(-)a + sin(-)0 , = - sin(-)a + cos(-)o , (29) 



For g = 1, one may resort to the random phase approximation (RPA) used 
in Ref.[33] to compute the "dressed" transverse spin susceptibility (plotted in 
FigJSD: 
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r«,Tr, r«»Tr -,(0) 

[X ] (£2) [X(£2)r 




r«,Tr -,(0) r^,Ti- -,(0) r„,Tr^ 

[X(n)r [Xio)] [X ] (H) 

Fig. 5. Diagrams in the RPA sum yielding Eq. ()30p . 




a 

Fig. 6. Plot of the spectral density S{Q,) in the various regime accessible to the 
device: Blue line: S{Q) near by the WFP; Red line: S{Cl) near by the SFP; 
Black line: S{n) near by the FFP (see Ref-IM] for details). 

where is given by 



[x^'f\Q) = j dr e^^^ (T.[{cos(e)S^(r) - sin(e)S^(r)}{cos(^)S^(0) - sin(0)S^(O)}])(o) 





4A 



^]2 + ■ 



(31) 



From Eqs. ( 12411311) ■ one gets 

'^m[x^'{-iVL + 0+)] 



S{9) 



4A 



4 



-^]2+4A2)2 



(32) 



A plot of S'(fi) vs. Vt is reported in Figj6l One sees that the limited broadening 
of the two peaks at f2 ~ 2A (which is the signal of frustration of decoherence) 
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depends on having both gi finite, and sin(^) 7^ 0. Furthermore, choosing g = I, 
imphes that gi is the coefficient of a marginal perturbation and, thus, is not 
renormahzed by the interaction. 



For g > 1, the system is attracted by an IR stable FFP and one has that gi 
goes to gl, where gl is the value of this coupling at the FFP. Even if in section 
[3] it has been shown that the FFP is IR attractive only when B\\ = B± = 0, 
the evaluation of the" dressed" transverse spin susceptibility requires to slightly 
move from the FFP by applying a small, nonzero B\\ and/or B±, in order to 
split the two impurity levels [33]. This may be safely carried out when the 
JJN has a finite size L since, in this case, the FFP is stable also against small 
fluctuations of the control parameters B±, provided that u/L is sufficiently 
big (see appendix [B|) . As a result, for g > I, the equation yielding S{Q) has 
the same form as Eq. fl321) . provided one replaces gi with gl and introduces 
a "self-energy correction" oc [Mj. From the duality relations given in 

Eq.f l22l) . one finds that Eq. (l32|) holds also for g = 4. Our results confirm that 
frustration of decoherence may be a remarkable signature of the emergence of 
a FFP in the phase diagram accessible to superconducting quantum devices 

EI 



4 Josephson current and coherent tunneling of Cooper pairs 



In section [21 we showed that the Hamiltonian i^B is invariant under the Ti- 
symmetry, even if it contains a term oc cos[$(0)]. In the following, we show 
that this symmetry is responsible for the appearance of 4e superconducting 
correlations in all the phases accessible to the JJN. For this purpose, we shall 
compute the Josephson current (JC) across C in all the phases of the JJN. 

To induce a JC, one may connect the leads of the JJN to two bulk supercon- 
ductors at fixed phase difference a (see Fig{T]) [10"13]. This amounts to require 
that, at the outer boundary (x = L), ^{L,t) = a. The (0-temperature) dc JC 
J(a) may then be computed as 



U„,£:^ . (33) 

/3^-oo fj da 



with e* = 2e being the Cooper pair charge and Z defined in Eq. 
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4-1 DC Josephson current near by the weakly coupled fixed point 



Near by the WFP, r) is given by 

a(n) \ f 1\ nx^ iWur 
$(a;,r) =a+^-^cos fn + -j— e-l"+2)— , (34) 



nez ^ + 2 



with 



Taking into account that (: e*"*'^°''^) :)(o) = e*"^", one gets 

In Z ^ In Zo + yi 1 rfr (: cos [$(0, r)] :)(o) + 1/2 / c/r (: cos [2$(0, r)] :)(o) 


= In Zq + cos(q;) + 1/2 cos(2a)} , (35) 

from which the dc JC is given by 

/(a) ^ e*yi sm{a) + 2e*y2 sin(2a) . (36) 



By inspection of Eq. (l36|) . one sees that /(a) is the sum of a term (oc sin(a)) 
corresponding to the usual tunnehng of Cooper pairs (of charge e*) across C, 
and of a term (oc sin(2a)) describing the coherent tunnehng of pairs of Cooper 
pairs (CTCP) (of total charge 2e*) across C. The relative weight of the two 
terms contributing to I (a) is oc |^ oc cos(6'). Thus, it may be tuned upon 
acting on till, at = vr, only the term oc sin(2a) is left, that is, the JC 
flows across C only because of CTCP. In FigJTl we plot the dc JC across C 
for different values of ip. 

The CTCP is also evidenced in the ac JC, arising across C when a finite 
voltage bias V is applied to the ends of both leads. To account for V, one has 
simply to add to i^BDSG a "voltage bias" term Hy = which yields a 

time-dependent shift of $ as $ — )■ $ + e*Vt, from which the ac JC is given by 

I{a, t) I{a + e*Vt) = e*yi sm{a + e*Vt) + 2e*y2 sin(2a + 2e*Vt) .(37) 
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Fig. 7. Plot of the dc Josephson current across C at the WFP for ^2/51 = -3 and 
for different values of B\\/B_\_. a): Top left panel: B\i^/B_\_ = 10; b): Bottom left 
panel: B^^/B± = .1; c): Bottom right panel: B^^/B± = .05; d): Top right panel: 
= 0. 
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^ 27U 7t 27t 



a 

Fig. 8. Plot of the dc Josephson current across C at the SFP for (from top left 
panel to top right panel, rotating counterclockwisely) Ae/|^ ^ 1; Ae/|^ ~ 1; 
Ae/fg « 1; Ae/fl = 0. 

4-2 DC Josephson current near by the strongly coupled fixed point 

Near by the SFP, the field given by Eq.f l38p . contains contributions from 
the zero mode operator. As a result, the partition function may be factorized 
as Z = Zo-modc["]Zosc5 where Zo_mode[tt] is the contribution of the zero-mode 
operators. From the analysis of the zero mode spectrum carried in subsection 
13. 2[ one sees that, for > 0, the eigenvalues of the zero- mode operator P are 
given by {pk} = {y/gik + a/vr)}, with the odd-k eigenvalues higher in energy 
by Ae = gi cos{9). As a result 
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Zo-mode[a] = Yl 

k&Z 



(3i 



from which the dc JC is obtained as 

e* (9 In Zo_modeN 



1(a) 



hm 



da 



(39) 



In FigJHlwe plot I{a) vs. a for different values of Ae/|^. One notices that for 
Ae/|^ ^ 1) the current takes the usual sawtooth behavior with periodicity 
in a equal to 2tt. Furthermore, as Ae/|^ ~ 1, one sees that the current 
continuously crosses a = n, with "satellite" jumps at a 



vr ± 



4LAe 



2 TTug 

As Ae/|^ approaches 0, the satellite jumps move to tt ± tt/2, where they 
eventually take place, when Ae = 0. As a result, the period of the sawtooth 
is exactly halved at Ae/|^ = 0. 



4-3 DC Josephson current near by the finite coupling fixed point 



To compute the JC near by the FFP, one must include the effects of the 
relevant SI perturbation at the SFP. While long instantons just provide a 
smoothing of the spikes of I (a) at the SFP [34J, for = and B± = 0, 
Sis become a relevant perturbation and drive the system to the IR attractive 
FFP. As a result, for = and B_i = 0, one expects that the net effect of 
the Sis on the dc JC will be to smoothen the sharp jumps of the sawtooth 
depicted in FiglHl while still preserving the half periodicity of /(a). 



In the following, we compute I{a) near by the jump at a = tt/2 for Ae/|^ = 
(top right panel of FigjH]). Upon tuning a so that < a < vr, the low energy 
states contributing to J (a) are associated with the zero modes labeled by 
Pk = Po and pk = —1. Their energies are respectively given by Eo{a) ± eo(a), 
with 



Eo{a) 



9 



AttL 



vr^ /vr \ 2 



eo(a 



9 



a 



(40) 



Since, for p = pq, = 1, while, p = = —1, one may trade the 

degeneracy-breaking energy eo(«) with an effective -B||(«) = eo(ct); in addition, 
as one is approaching the FFP from the SFP, one should take, as boundary 
interaction, the dual boundary Hamiltonian given in Eq.f l2ip : the term 
(oc cos[6(0)/2]) now acts on the spin S as an effective field B± ~ /2i. As a 
result, one may compute I{a) for a ~ 7r/2 by taking the logarithmic derivative 
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Fig. 9. Sketch of I (a) vs. a near by the FFP (Eq.([43])). 
with respect to a of the relevant contribution to Zo-mode[tt]; namely 



(41) 



with 



SE[{a,b}] = J dr 



d_ 



B\\{a) 



d_ 



Bu(a) 



while the ellipses denote contributions coming from pk 7^ {PoiP-i}- Using 
Eq.f l33|) . one gets 



Ka 




(a-f) +(^) L-f 



(43) 



Eq.f l43|) shows that /(a) is a smooth function, linearly crossing 0, at a = 7i/2. 
The term oc L^~2 in the denominator evidences that this behavior is more 
pronounced, as the system size increases. Carrying out a similar computation 
for a 7^ 7r/2 leads to the plot displayed in FiglQl 

In the next section, we shall show how the results obtained for the coherent 
current transport in the JJN extend to the situation in which a dissipative 
current is induced, as well. 
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5 DC transport and coherent tunneling of Cooper pairs 

In this section, we analyze the observable consequences of CTCP on a dc trans- 
port experiment. For this purpose, we shall compute the dissipative current 
/dc(V^) flowing across C when the leads are connected to two metallic wires, 
at flxed voltage bias V. Differently from the Josephson supercurrent, /dc(^) 
fluctuates even at T = 0, in presence of weak tunneling and/or backscattering 
centers, and this zero-temperature fluctuations are measured by the shot noise 
at voltage bias V and frequency fl, 



where {, } is the anticommutator, (. . .)ne denotes averaging, at zero tempera- 
ture and at flnite V, la is either identifled with the transmitted current, when 
the two leads are weakly coupled to each other, or with the backscattered 
current, when the leads are strongly coupled. 

Measuring the shot noise provides a tool to probe the effective elementary 
charge q flowing across a device since, from the functional dependence of 
S{V, Q) upon V, one obtains q as 



Measurements of S have been used to probe noninteger charges in fractional 
quantum Hall bars [28], in normal metal-superconductor junctions [29], in 
Kondo dots [27J. In the following, we show that, by properly tuning the exter- 
nal control parameter, pairing of Cooper pairs may be as well detected in a 
shot noise measurement, provided that V is lower than a threshold, depending 
on the values of B\\ and B±. 

In subsection 15.11 we compute I&iV) and S{V) at flnite V near by the WFP. 
We show that, when B\\ =0, B± deflnes a threshold below which the effective 
charge detected from a shot noise measurement is g = 2e*, which evidences a 
phase with 4e superconductivity. 

In subsection 15.21 we use the fermionization procedure of the BDSG Hamilto- 
nian ioi g = 2 [36|, to extend the results of subsection 15. II beyond the pertur- 
bative regime in the boundary couplings. Mathematical details concerning the 
derivation of subsection 15.11 are provided in appendix [Cl while the fermioniza- 
tion procedure is discussed in detail in appendix [Dl 




(44) 



q{V) = S{V,n^0)/i2hiV)) 



(45) 
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5. 1 Perturbative computation of the dissipative current and of the shot noise 
at finite voltage bias V 

When the JJN is connected to two metaUic wires at finite voltage bias V, 
the computation of the current flowing across C -and its fluctuations- may be 
carried out within the Hamihonian formahsm. In particular, near by the WFP, 
one may consistently compute both /dc(V") and S{V) nonperturbatively in V, 
resorting to time-dependent perturbation theory in the boundary couplings gi, 



Applying a finite voltage bias V amounts to shift the field $ as $ — )■ $ + uyt, 
with uv — e*V. To the leading order in 5^1,5^2, the current operator in the 
Heisenberg representation, jnit), is obtained by 



where ji{t) and H^{t) are, respectively, the current operator and the boundary 
Hamiltonian in the interaction representation, given by 

ji{t)^e*giS^{t) : sm[^{t) + uyt] : +2e*g2 : sin[2$(i) + 2ujvt] : 
Hs(t) = -giS^{t) : cos[$(i) + uyt] : -^2 : cos[2$(f) + 2ujvt\ : , (47) 

with $(t) = ^{x = 0,t). 

The 0-temperature average value of jnit) is given by 




(46) 



—00 






(48) 



while S{V) is given by 




23 



-2(e 



*\2 • 2 

^ sm 



dr 



r TTiU 



e i 



/ rf^ 1_ 



(49) 



The relevant integrals needed in Eqs. fl48ll49p are computed in appendix O As 
a result, /dc(^) and SiV) are given by 



/Dc(V^) = e* cos^(^^)^f3fte z 



r[i-!]r[i^ + i + i] , r[i-f]r[-^ + i + i] 



+ 



K-^)r[^-i + i] (u;y + ^)r[ 



2ttu 



+e sm 



in2(^) ^eU 



r[i - f] n-^^^ + i + 1] 



1] 



r[i-f]r[^ + f + i] ^ r[i-^]r[- 



+ ^11 



{2u^ _ |m) r[iHf - 1 + 11 (2..,, + |y) r[- 



1] 



(50) 



and 



5(\/) = 2(e*)"cos'(^)^f g^e 2 



r[i-f]r[-^ + i + i] 



+2(e*)^^f sin 



3in2(^) sRej? 



r[i-f]r[ 



L(a;v+2A) 
2ttu 



+ 1] 



2nu 



1] 



27rM 



[uv - 2A + ^) rh ^^"'-'^^) 

+ 8(e*)2^^ 3?e(z 



1] 



+ 



r[i-^]r[- 



Lu)v 



+ 1 + 1] 



- f + 11 + ^) r[- 



1] 



.(51) 



Finally, the Stirling's approximation (Eq. flC.5P ) allows to derive the large-L 
limit {{Luv/2nu) > 1) of Eqs. fl50]l5T|) . leading to 

I^ciV) ^ e r[l - -]sm ( - n [-^) 



9 



9 
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+ 



Uy — 2A 



9{uv - 2A) 



' L{tuv-2\) y 



g \ g I 2uv \ tcu 



8 

V ^ » 



(52) 



and 



gfsm\9) 



+ 



OJy — 2\ 



9{uv - 2A) 



g \ g J \ OJy V 271U J 

'L{uJv - 2A)\ ' 



2'KU 



+ 8(e*)'r[l - -] sin — 



7r\ gi [Lux 



g J 2uv \ vr-u 



,(53) 



with 9{x) being Heavside's 6'-function. 

To determine qiV), one observes that the term cx ^fsin^(6') in Eqs.( l52|53|) 
exhibits a threshold at cuy = 2A. In computing qiV), one notices that, since 
the term oc cos[<l'(0)] and the term oc cos[2$(0)] in have different scahng 
dimensions {1/g and i/g, respectively), tuning the low-energy scale uy induces 
a flow in q{V). Indeed, assuming B\\ ^ 0, one finds that qiV) ~ e* for <^ w*, 
while qiV) ~ 2e* for uy ^ u*, with the crossover scale given by 



id* 



T 



glcos\9)ni 



9l 



r[i 



sm 



sm 



(f) 



(54) 



Eg. (15^ states that the contributions to /dc(^) and to SiV), which are oc 
gfcos'^{9) and oc (?| respectively, are of the same order of magnitude. 

As Uy — )■ 0, the value of q{V) characterizing the IR stable fixed point is always 
q = e*, except when B\\ = 0,Bj_ ^ 0, where, the existence of a threshold at 
Uy = 2A implies that ^(V") = 2e*, for uy < 2\ = 2|i?^|. Thus, for _B|| = 
0, B± 7^ 4e superconductivity is evidenced in the dissipative dc transport as 
well, when uy < 2A. qiV) is plotted in Fig{TUl both for 7^ 0, and for _B|j = 
0, B± 7^ 0. For B\\ = 0, B± = 0, the threshold at uy = 2A moves to uy = 0. 
This is again consistent with the picture of a 4e-deconfined superconductor, 
with 2e elementary charged excitations (see section [3?2|) . 



5.2 Exact formula for /dc(^) o,nd S(y) in the fermionized theory for g = 2 



In this subsection, we use the exact fermionization of the BDSG Hamiltonian 
at g = 2 (see appendix [D] for details), to compute /dc(^) and S(y) beyond 
perturbation theory, for B» =0. 
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Fig. 10. Plot of q{V) = S{V)/Idc{V). Left-hand panel: q{V) vs. ioy for B\\ / 
(ri -symmetry explicitly broken); Right-hand panel: qiV) vs. ujy for B\\ = 
(ri-symmetric case). While for B\\ ^ q{V) flows from e* (low- voltage limit) to 
2e*, for B\\ = and S_l / q{V) = 2e* also for |a;y| < 2A = 2|S_l|. 

Applying a finite voltage bias V to the leads amounts to add to the fermionic 
action S^^^ given in Eq.( ]D.8l) . a contribution given by 



Sv = ^ [ dt [ dx[■.^P^Ji{x,t)Mx,t)■ + ■■^L{x,t)M^,i)■] ■ (55) 







In the fermionized theory, the current operator jnit) is given by 

jH{t) = e*u{: ^Jj(0, t)^K(0, t) : - : ^^(0, t)MO, t) :} . (56) 



Taking into account that the finite voltage bias shifts the chemical potential 
of the R ( l) modes by e*V/2 {—e*V/2) and using the linear relations between 
the right-handed modes iPr{p) and the left-handed modes ipiip) provided in 
Eq. (]D.ll[) . one obtains that the average current /dc(^) is given by 

{e*V/2u) 

lT,dV) = — [1-\R{P)? + \Q{P)?] = / dp\Q{p)\' ,(57) 

^ 0<p<{e*V/2u) ^ 



where, to derive Eq.(!57l). we use |Q(p)P = 1 — |-R(p)p. 

To compute SiV) through Eq. fH^ . one first needs to determine S(t,t') = 
({/fl(t), /fl(t')})NE, which, in the fermionized theory, is given by S(t — t') = 
Si(t - t') + Ssjt - t') + Esit - f), with 

Si(t - f) = -2 r^)' 5]e^[^.-.K*-0[i _ R*(p)R(q) + Q*(_p)g(_g)] x 
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(58) 



E2(t-0 = S3(t-t') 



and ep = up. 

Taking the Fourier transform of S(t — t') at frequency Q and setting = 0, 
one obtains 

S{V) = S{Q^O,V) = ^-^ Y: \Q{p)\' \R{p)\' = 

-{e''V/2u)<p<{e*V/2u) 

{e*V/2u) 



]ie*yu 



J dp {\Q{p)\' - \Q{p)\'} . (60) 



As a result, from Eq. (lD.13P of appendix [Dl one has 



(e*V/2u) ^ „2 „ _ 

T n/^ 2e*M f [-^ + ug2P^ + 2B^g2pf 
Idc{V) = / dp ^ -2 , (61) 







and 



SiV) = ?<£:£^ "T' dp (-V+^B,n-i + -9-^^^^B,gM' .(82) 



Both integrals in Eqs. fl61f62p may be exactly computed, although the final 
expressions are not, in general, very enlightening. However, for = 0, one 
may plot q{V) vs. V, obtaining a result matching the left-hand panel of FigJTU] 
in particular, one finds that 

lim q(V) = e* , lim q(V) = 2e* . (63) 



Eq. (l63!) shows that the crossover between an effective elementary charge e* 
and an effective charge 2e*, induced by acting on the applied voltage bias V, 
is generally valid, independently of the reliability of the perturbative compu- 
tation of subsection I5.1[ 
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For B± 7^ 0, one obtains 

lunq{V) = , 32e*(2S^)^ ^^^^ 

16 (li) + ^ (f ) (2|5J) + (16 + (f ) )(2i?^)^ 



which, to the fourth order in gi,g2, yields 

hm q(V) = 2e* . (65) 
y-5-o 

This is again consistent with the resuh of Eq. fl53|) . since it shows that B± ^ 
may stabihze an IR stable fixed point exhibiting 4e superconductivity. 



6 Concluding remarks 



We investigated a boundary double Sine-Gordon model, where two boundary 
terms, of different periodicity and scale dimensions, are coupled to a Kondo- 
like spin degree of freedom. We showed that the pertinent engineering of the 
coupling between the spin degree of freedom and the bosonic field induces 
a competition between the two boundary terms, and that this gives rise to 
nonperturbative phenomena, such as the emergence of novel quantum phases. 
We showed indeed that the strongly coupled fixed point- in a pertinent range of 
parameters- becomes unstable as a result of the deconfinement of new phase- 
slip operators (i.e., the short instantons), arising from adding the less relevant 
boundary operator. 

To look for a physical context where such nonperturbative phenomena may be 
observed, we analyzed a Josephson junction network, providing a remarkable 
realization of the BDSG field theory described by Hb- For this network, we 
showed that the competition between the two periodicities stabilizes a robust 
finite coupling fixed point and allows for the emergence, at a pertinent scale, 
of 4e superconductivity. To probe the latter phenomenon, we computed the 
dc shot noise and showed that, by pertinently tuning the applied voltage bias 
V, the effective charge of the carriers varies from 2e to 4e. 

In our analysis the onset of the ti symmetry is tuned by acting upon the 
control parameter . The role of this symmetry is twofold. In section 13.11 we 
showed that, for B± 7^ 0, 4e superconductivity is naturally associated to the 
realization of this symmetry. On the other hand, in section [3l2] we showed that, 
for B± = 0, Ti-symmetry implies SI deconfinement (and, thus, for 1 < g < 4 
it stabilizes the IR stable FFP); at variance, short instantons are confined, if 
Ti-symmetry is "broken", and the SFP is IR stable for 1 < g. 



28 



We argued that while, as the system size L — )■ oo, the FFP is IR stable only 
when B\\ = B± = 0, for a realistic JJN of finite size L, the FFP is robust 
also against small fluctuations of the control parameters B\\,B_i_, provided 
that u/L is sufficiently big. Indeed, as discussed in detail in appendix |Bl SI 
deconfinement (which is the mechanism de-stabilizing the SFP vs. the FFP) 
is effective only for L > Lstop ~ cos(6')). As a result, for L < Lstop, 

one may safely assume that the behavior of the device is still driven by the 
FFP. This renders the emergence of the FFP in the phase diagram relevant 
for engineering realistic superconducting devices with enhanced coherence. 

It may be worth to investigate how the methods developed in this paper could 
be modified to account for the effects of different commensurability ratios 
between competing boundary interactions of a boundary Hamiltonian. 
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A Derivation of 



In this appendix, we describe the Schrieffer- Wolff summation procedure, to 
derive the boundary interaction Hamiltonian in Eq. ([2]) . In order to imple- 
ment this procedure, one should determine the low-energy states of the central 
region, obtained after diagonalizing the effective spin-1/2 Hamiltonian 

Hc = -jj2{e'^S;s]^, + h.c.}-hj2S^ . (A.l) 

j=0 j=0 

In Eq.f lA.l|) . h describes a slight detuning of the gate voltage Vg = N + 1/2 
{N integer), acting on the junctions in C. 



A.l Eigenvalues and eigenstates of Hq 



He commutes with = Z]j=o ^jj thus, its eigenstates may be grouped in 
multiplets of 5*^. Namely, the spectrum is given by: 

• Spin-2 eigenstate There is only one state, of energy 62, given by 

|2) = |tttt) , €2 = -4/1 , (A.2) 

• Spin- -2 eigenstate There is only one state, as well, of energy e_2, given 
by 

|-2) = |iiii) , e_2 = 4/i , (A.3) 

• Spin-1 eigenstates There are four spin-1 eigenstates, with energies given 

by 

\l,k) = lj2e'''\l,j) , e^{k) = -2Jcos{k-j)~2h , (A.4) 
where 

27r£ , , 

k=— , A; = 0,l,2,3 . (A.5) 

is the spin-1 state where all the spins are except the one at site j, 
which is l- 

• Spin- -1 eigenstates There are four spin-(-l) eigenstates as well, with 
energies given by 

|-l,fc) = if]e*'=Vl,j) , e.,{k) = -2Jcos{k + j) + 2h , (A.6) 
^ i=o ^ 
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where 



, /c = 0, 1,2,3 . (A.7) 

I — 1, j) is the spin-1 state, where all the spins are J,, except the one at site 
J, which is 

• Spin-0 eigenstates There arc six spin-0 eigenstates. They are listed below, 
with their corresponding energies 
State |0, a) 

|o, a) = ^{V2[| tm> + 1 mt)] - 1 utt) 
-im;)-imt)-iim)} , e, = 2v^jcos(|) , (a.s) 

state 1 0,6) 

|o, b) = ^{V2[\ nu) - 1 mt)] + ^l un) 

+ i\nU)-i\nW-i\mi)} , e;. = 2^27 sin , (A.9) 

state |0,c) 

|o, c) = ^{ V2[| nu) + 1 mt)] + 1 un) 

+ im;) + imt) + um)} , ee = -2^jcos(|) , (a.io) 

state \0,d) 

|o, d) = ^{V2[\ nu) - 1 mt)] - i\ un) 

-i\nU)+i\nW+i\mi)} , e, = -2v^Jsin(^) , (A.ll) 



state |0,e) 



o,e) = i={im;)-i;;n)} , e, = o , (A.12) 



V2 

state |0,/) 



|0,/) = -^{|mt)-Um)} , ef^O . (A.13) 
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For = TT, the groundstate is twofold degenerate: both states \0,c),\0,d) 
have the minimum possible energy (—2 J). In the next section, we consider 
the effective boundary Hamiltonian arising when the states are coupled to the 
leads. 



A. 2 Coupling C to the leads: effective boundary Hamiltonian 



To determine the effective boundary Hamiltonian describing the low- 
energy dynamics of C connected to the leads, one starts by connecting C 
to the leads by two junctions, of nominal strength A <^ Ej, J, with the Hamil- 
tonian 

i^T = -A{e^*<^°^^o- + e^*>^°^^2" + h.c.} , (A.14) 



where So, S2 have been defined in FigJH For y9 ~ vr, one should consider the 
low-energy doublet spanned by the states |\l'i-) = |0,c) , I'^i) = \0,d), and, 
in order to implement the SW summation, one needs the matrix elements of 
|\E'-f-) or |\E'j^) with the lowest-lying excited states of He- The relevant matrix 
elements are listed below: 

• First set 

(l,0|ifT|^4) = ^[e-75*<('') +e-75*>(0)] ^ 

• Second set 

{-1,0\Ht\^^) = -^[e73*<(0) +e7f*>(0)] ^ ^^^^g^ 

• Third set 

(1, l\H.\^,) = -l^±^[e-7I*<(0) - e-7f*>(0)] , (A.17) 
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• Fourth set 

{-i,||ffT|*t) = 3l^**^'°'-^**^'°'l 

(-l,||ff,|*,)=<^l^|e*»<'»'-e**>'»'l . (A.18) 

Using the matrix elements given in Eqs.( ]A.15|IA.16|IA.17|IA.18p . the SW sum- 
mation yields, to second order in A, the following contribution to the effective 
boundary Hamiltonian 

\2 r r ^*<(o) , ^'i'>(o)ir -^'i'<(o) , — 7=*>(o)i 

Etl = |vl/t)(^tl-|(2 + V2f ^' \ ^ 

Eff I t/\ tl32\v J _(2_72)J + 2/i 



+ (2 + V2 



+2- 



+2. 



-(2 - y^)J-2h 
-(2 - V2)J- 2/i 



32 r ' -(2-V2)J + 2/i 

+ (2 + ^2)^^ ^ 



+2 



+2 



-(2 - V2)J -2h 
-(2 - v^)J + 2/i 



-[|vl/^)(vl/j + |v[/^)(vl/^|] 



(2- 72)7-2/1 

v/2(2 + v^)A2 r [e75*<W + e7f*>(°)][e-7f*<W + e-7f*>(°)] 



32 I _(2-V2)J + 2/i 



-(2-^2)7-2/1 
[e75*<(0) _ ^-^<i.,(o)j^^-^<i.<(o) _ g-;^<i>>(o)j 



-(2 - ^/2)J + 2h 
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-(2 - V2)J + 2/i 



(A.19) 



Notice that, in Eg. flA.lQp . = tt. A slight detuning of off = tt just adds 
to Hj^g and effective energy sphtting term, given by 



SHii = 4 J sin ( ) |^t)(^J = -^11 



(A.20) 



where, in Eq. (IA.20p . one sets 5|| = -4Jsin(^^) and S"^ = - 
Expanding Eq. (lA.19P to the leading order in the detuning parameter 
h, one finds 



2J \ 2-V2, 



$>(0)-$<(0) 



V2 



(2 - V2)2 7^' 



(A.2f) 



with = |\E'^)(\I^j^| + The BDSG Hamiltonian is recovered by 

following the same procedure, up to fourth-order in A. The additional term is 
given by 

H^) = £^[e75[*<W-*>Wl - e-7l[*<(°)-*>(°)l]2 ^ ^ (^.22) 



where C ~ 10 ^ is a numerical coefficient and the ellipses stand for subleading 
contributions. Defining $ = ($<- ^y)/y/2, gi = ^ (2^)' ^2 = 2^ and 

B± = ^^2--/Ty^^ '^^^ obtains the effective boundary Hamiltonian given in 

Eq.(|2]), as i^B = -f^Eff + -^Eff + ^-^11- We notice that both B± and i^y may be 
tuned by acting on external control fields: B± with Vg, B\\ with if. 



B Instanton solutions of the boundary double Sine-Gordon model 
and short instanton deconfinement 



In this appendix, we study in detail instanton solutions of the BDSG model, 
regarded as imaginary time trajectories of the zero mode, -P(r). To construct 
the effective Euclidean action for P(r), Se[P], one uses for r) the mode 
expansion in Eq. (fT9|) . S'£;[P] is computed from 

e-s^[P] = J D$„,,D{a, b} e'^'l^l-^B ^ (B.i) 
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where / D$oscD{o, b} denotes functional integration over the oscillator modes 
of the field $, as well as over the local fermionic operators, associated with 
the spin S. The free Euclidean action for the field r) is given by 



Sf[^ = 1- J dr J dx 



1 /d^' 
u \ dr 



+ u(— 

\ dx , 



(B.2) 



From the the mode expansion given in Eq.( ]T9|) . one obtains 
SE[P] = sf[P] + 5SE[P] , 



(B.3) 



with 



4°' IP] 



n 

uL 



— gi cos(0) cos 



TT 



-P -a 



92 cos 



27r 

71^ 



-p 



M = nL/6u, and 6Se[P] defined by 



(B.5) 



From Eq.( lB.4l) . one sees that the inductance energy (oc L"^) breaks, in the 
finite-size system, the degeneracy between the minima of H^] however, a de- 
generacy between only nearest neighboring minima of may be restored by 
setting a = TT + 2kn, if B\\ > 0, or a = 27rfc, if B\\ < 0. 

Assuming > and a = — tt, a single-instanton solution, representing a 

i55'(o) rpi 

quantum jump between pq and p_i may be built by requiring that — fp-^ = 0. 
Defining p(r) = (-P(t) — l/2)/y^, one obtains the following imaginary time 
"equation of motion" for p{t): 

Mp ^{p — tt) + T^gi cos{9) cosfvrp] + 27ig2 cospvrp] = . (B.6) 



Apart from the term cx Eq. flB.6P is the imaginary-time version of the equa- 
tion yielding static, finite energy, soliton solutions in the double Sine-Gordon 
model |j32j. Borrowing well-known results [32\, one may write down a single- 
instanton solution as 



p{t) = — 1 H — -i arctan 

TT 



exp 



T + R 



— arctan 



exp 



-T + R 



.(B.7) 
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In Eq. (lB.7p . the "bare" parameter R is defined by the condition ^ sinh^ (jji cos{9) + R 
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gicos{ey analysis, due to the logarithmic divergences induced by inter- 

action with the oscillator modes, a logarithmic term in R is induced. Due to 
this, we regard i? as a variational parameter and use the "instanton size" T, 
as the main scaling parameter. 



From Eq. flB.7p . p{t) may be represented as a sequence of two short instantons, 
separated by a distance 2R from each other. In addition, Eq. flB.7p allows to 
compute also the Euclidean action for p{t), yielding 



So^ -g\n{ — 



fuT\ , /uR\ AgR 



(B.8) 



from which, one gets 



e ^° (X. R^ exp 



AgR 

TTUL 



2gi cos{e)R 



(B.9) 



Eq. (IB.9p implies that large i?-solutions may either be favored, or disfavored, 
according to whether > 1, or < 1. Since ~ L'^~^, for g > 1, 

one sees that solutions with the two Sis at large separations are strongly 
suppressed, as L goes large. The optimal value of the variational parameter 
R is set by requiring that ^ = 0, which, for for small ^icos(6'), leads to 
-R* oc 1/gi cos{9). 

One should notice that, since R is directly related to cos(6')5fi, one may readily 
tune it by just acting on the applied fiux (p. Since Sis are confined over a 
scale Lstop ~ R*, the scaling of the parameter Ai(L) in Eq. fl23l) will stop at 
L = Lstop- To show how acting on ip may trigger SI dconfinement, in Fig JB.ll 
we plot the solution p(r) for different values of ip. The SI deconfinement may 
clearly be seen, for ip = tt, 

Finally, we mention that the solution for B\\ < and a = 0, has the same 
form as the one in Eq. (lB.7p . provided that one substitutes B\\ with 



C Tables of relevant integrals 



In this appendix, we sketch the calculation of the integrals used in subsection 
IS.ll to compute /dc(^) and S{V) near by the WFP. All the relevant integrals 



36 



P(x) 





(p=l.l7i; 



(p=1.0l7l 





X 



Fig. B.l. Instanton profile p{t) for various values of (p. From top to bottom: 
93 = 27r, l.lvr, l.Olvr, vr. The plots have been vertically shifted: in any case p{t) 
ranges from to 1. The dashed portion of the plot for if = tt (black line) denotes 
infinite separation between the short instantons (i? — 00). 

may be recasted in the form 



with ( real and rj = 0^. To compute the integral, one assumes z — Qt^u/ L < 0, 
and integrates over the closed integration path shown in Fig JC.ll to get 



00 




— IZT 



(C.l) 




(C.2) 



The integral in Eq. (lC.2p is tabulated [38], yielding 




(ttu 



L 



(C.3) 
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Im(x) 



Re(T) 




Fig. C.l. Integration path used to compute I(^{z) in Eq. llC.lj) 



where 2F1 is the confluent hypergeometric function. Then, using the identity 



in Eq. (]C.3p . one obtains Eqs.( [50|5T|) of subsection 15.11 To analytically extend 
Eq. (lC.3l) . one may use Stirling's approximation, which yields 



Eq.( ]C.5|) has been applied in subsection 15. when extracting the large-L limit 



D Fermionization and exact solution of the BDSG model for g = 2 

In this appendix we carry out the fermionization procedure of Hbdsg = 
HhL + Hb for (yf = 2. As we will show, introducing a complete set of fermionic 
coordinates to represent the relevant fields of the model, allows for recasting 
Hbdsg in ci quadratic form, thus making it exactly solvable. In the following, 
we assume = and, due to our interest in nonequilibrium dc transport 
properties, we resort to the real time formalism. 

The fermionization of Hbbsg follows the same basic steps as in Ref.[36j. To 
begin with, one fermionizes i^LL for = 2 by writing $ as the sum of two chiral 



2F1 [a, b, c; 1] 



T[c]T[c-a-b] 
T[c-a]T[c-b] 



(C.4) 




(C.5) 



of /dc(^) and ofS{V). 
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fields, (f)R,(f)L, as = (t)R{x — ut) + (t)L{x + ut). The chiral vertex operators 

■ ^-i(t>R{x-ut) . ^-i<j>L{x+ut) . j-Qg^y ^j^g regarded as two chiral fermionic fields 



iPr{x - Ut) = T]R 



1 



: e 



L 



where rjR, rji are the Klein factors, introduced to ensure the correct anticom- 
mutation relations between i/jr and ipL 



To fermionize H^, one assumes Neumann BCs at both boundaries ^6j; namely 
MO,t) = MO,t) , ML,t)=ML,t) . (D.2) 



Eg. (ID. 2 p is enforced if one uses [36] the noninteracting real time action for 



oFer 
'-'0 



i dt dx \ ipRiO,t) 



d_ d_ 

dt dx 



^fi(o,t) + v^i(o,t) 



d 



d 



dt ^ dx 



^L(0,t) 



+ i| / dt {4(0,t)^i(0,t)-^I(0,t)7/;^(0,t)} 



(D.3) 



Then, one has to follow a different procedure for the term oc gi and the term 
oc g2 in H^. Fermionizing the former, requires introducing additional local 
degrees of freedom describing the spin- 1/2 variable S emerging at C. 

Using a "rotated" Jordan- Wigner (JW) transformation [31], one introduces a 
complex fermion ^, in terms of which S^, are given by 



(D.4) 



Using Eqs. (ID.l"|D.4p . and taking into account the boundary conditions in 
Eq. (ID.2p . one gets the contribution to boundary action which is oc (71; namely 
0. 



rft{[4(0,t)+^i(0,t)]e-e^[^i^(0,t)+V^I(0,t)]} + 2 / dtS}^-^ .(D.5) 



dt 



To fermionize the contribution to the boundary interaction which is propor- 
tional to (72, one has to regularize, by point splitting, the products : e"*'^'*^^^ : 



Notice that gi, not gi, appears in Eq. ()D.5p . as the scale dependent factor has 
been reabsorbed in the definition of the fermionic fields, Eq. ()D.2p 
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■ g6i<AB(o) -with a,b = ±1 and A,B = L, R, and to require that the anticom- 
mutation relations between the chiral fermions are preserved. As a resuh, one 
gets 



^(11) _ ^92 



dt <^|j(0,t) — + V^L(0,t)- 



dt 



dt 
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c/t <;v^)j(o,t) — - — + ^h{o,t) — - — 



dt 



(D.6) 



Finally, the term oc Bj_ is fermionized using the JW transformation reported in 
Eq.( 1D.4p and, by adding the "kinetic" term for ^, one gets the last contribution 
to the boundary action, which is given by 



qini) 



dt 



1 



(D.7) 



From Eqs. (ID.3|D.5|D.6|D.7p . one sees that the full fermionic action is given 

by 



r.Fer _ r.Fer q{I) . q{II) . q{III) 



(D.8) 



Equating to zero the functional derivative of S'^'^'', one obtains the boundary 
conditions for the fermionic fields, which are given by 

(f^ +\ I tn , -^2 f d'ipl{0,t) d^R{0,t) \ gi 

2 [^^(0' - MO, t)] + r- [—Q^ + j + ^y^ = 

^e + 5±e-^|[^/^(0,^) + ^i(0,t)] = . (D.9) 



Getting rid of in Eqs.( ]D.9"|) and using the normal mode expansion for 
i^R,ipL, given by 

V L p y/ L p 

one eventually obtains the following linear relations between the normal modes 

i^R{p),i^L{p) 

= R{p)i^L{-p) + QipWhip) 
^li(-p) = Q{p)M-p) + R{p)A{p) , (D-ii) 
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with 



up + 2B\ ,^ 
R{p) = -il^ ^- , (D.12) 



and 



up + 2B^-i[-^ + ug2P^ + 2B^g2p] 



Eqs. (1D.12|D.13P provide the formulas we used in subsection 15.21 to compute 
Jdc(^) and S{V) in the fermionized theory. 
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